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Multiwavelength afterglow light curves from magnetized GRB flows 
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ABSTRACT 

We use high-resolution relativistic MHD simulations coupled with a radiative transfer code 
to compute multiwavelength afterglow light curves of magnetized ejecta of gamma-ray bursts 
interacting with a uniform circumburst medium. The aim of our study is to determine how the 
magnetization of the ejecta at large distance from the central engine influences the afterglow 
emission, and to assess whether observations can be reliably used to infer the strength of the 
magnetic field. We find that, for typical parameters of the ejecta, the emission from the reverse 
shock peaks for magnetization cr ~ 0.01 - 0.1 of the flow, and that it is greatly suppressed for 
higher <xo. The emission from the forward shock shows an achromatic break shortly after the 
end of the burst marking the onset of the self-similar evolution of the blast wave. Fitting the 
early afterglow of GRB 990123 and 090102 with our numerical models we infer respective 
magnetizations of cr ~ 0.01 and ctq ~ 0.1 for these bursts. We argue that the lack of observed 
reverse shock emission from the majority of the bursts can be understood if <x > 0. 1 , since we 
obtain that the luminosity of the reverse shock decreases significantly for ctq ~ 1. For ejecta 
with cto > 0.1 our models predict that there is sufficient energy left in the magnetic field, at 
least during an interval of ~ 10 times the burst duration, to produce a substantial emission 
if the magnetic energy can be dissipated (for instance, due to resistive effects) and radiated 
away. 
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1 INTRODUCTION 

Gamma ray bursts (GRBs) are thought to be caused by the en- 
ergy dissipated by an ultra-relativistic outflow. Two alternative 
scenarios for powering the flow have been extensively consid- 
ered, resulting either in the ther mal-energy dominated fireball 



jGoodman|[T98r3 : 



Paczynski 1986) or in the Poynting-flux domi 
Usovlll992l ; lThompsorj[l994l : iMeszaros & Reesl 



nated flow (PDF: _ 

1 1997b . In the fireball model, magnetic fields are not dynamically 
important, and at large distances the flow is weakly magnetized 
<x <k 1. Radiation pressure drives the acceleration of the flow 
until the bulk Lorentz factor becomes To - h m (where h m is the 
initial specific enthalphy of the flow; see e.g.. lAlov et alj|2005l : 
I Aloy & Obergaulingeil |2007|). and intern al collisions power the 
prompt emission jRees & Meszaroslll994l) . MHD models for GRB 
jets assume that the flow is launched Poynting-flux dominated 
(tr in » 1 ; just as in the widely accepted theoretical model for rel- 
ativistic jets from active galactic nuclei). The acceleration of the 
flow is model dependent and therefore uncertain. PDF flows have 
in common the fact that the acceleration to ultra-relativistic speed 
happens b y converting pa rt of the magnetic energy into kinetic en- 
ergy (e.g.. lLi et al. |[l992|). However, there is no full consensus on 
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which is the ma gnetization of the flow (ctq) when the accelera- 
tion is completed. iLvutikov & Blandfordl ( l2003h consider a sub-fast 
magnetosonic flow that maintains tr » 1 up to distances of r > 
10 16 cm, where it interacts with the surrounding medium. In con- 
trast, a large number of ideal and non-ideal MHD models for jet ac- 
celeration are super-fast magnetosonic. Still, they support a picture 
where the acceleration/dissipation processes are not 100% efficient 
in consuming the magnetic energy and, therefore, the flow magne- 
tization <t ~ 1 and t he bulk Lorentz factor is Tp ^ Q" in /(1 + q- ) ~ 
ir in at large distances jDrenkhahn & Spruitl2 002: Komissar ov et al.l 



120091: iTchekhovskov et al.ll2009l ; lLvubarsky||201Cl). During or af- 
ter the acceleration process magnet ic dissipation iThompsonll 19941: 
Spruit et al.l l200ll; iGiannios! 20081) or internal shocks dFan et al.l 



2004bl : iMimica & Alovl 120101) may power the GRB. Here we fo- 



cus on the dynamics of the deceleration of cr < 1 flows with focus 
on the role of magnetic field strength to the resulting emission. 

If the magnetic fields of the ejecta are dynamically unim- 
portant (e.g., (To <k 1) at large distance where the interaction 
with the circumburst medium becomes substantial, a forward shock 
(FS) and a reverse shock (RS) are expected to form at the inter- 
face between the G RB ejecta and circumburst (external) medium 
dSari & Piranlll995h . If the ejecta is strongly m agnetized, we ex- 
pect a weak or absent RS ( IGiannios et al. H2OO8I) . Observationally, 
there is a wealth of recent evidence that magnetic fields shall be 
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dynamically important. Most bursts show no trace of the RS emis- 
sion, in contrast to predictions of the fireba ll model (see, however, 
INakar & PirarJ 20051 : IMcMahon et alj2006l) . When the bump in the 
early optical light curve is observed (and resembles that expected 
by the RS), modeling of the light curves indicates that the magnetic 
field is l arger in the shocked ejecta than in the shocked external 
mediu m jFan et all 1 20021 : Izhang et al. 20031: IPanaitescu & Kumar! 
|2004) . More conclusive is the detection of ~ 10% polarization 
in the early optical afterglow emission of GRB 090102, consis- 
tent wi th the ejecta contai ning quasi-coherent, large-scale magnetic 
fields jSteeleetal.1 [20091) . Since dynamically dominant magnetic 
fields would suppress the RS emission altogether, these authors 
give a rough estimate of cro ~ 1 for GRB 090102, supporting super- 
fast magnetosonic MHD models for GRBs. 



To make quantified inferences for the degree of the magnetiza- 
tion of the ejecta, a detailed understanding of the dynamical inter- 
action between magnetized ejecta and external medium is needed. 
This has been done at an increasing level of sophistication. As a 
first approach to the problem, the dynamical effects of magnetic 
fields have been ign ored, assuming a purely hydrodynamic be- 
havior of the ejecta teees & Meszaros! Il992l; ISari & Piranlll995l ; 
iKobavashi etail 1 19991 : iGranot et al.l l200ll : iMeliani et alj 120071) . 
Such an approximation has lead to infer that the magnetic field 
strength at the RS and at the FS is different (i.e., the parameter 
e B is larger at the RS than at the FS jFan et aljr2002l : IZhang et al.l 
2003)). A more accurate approach takes into account the mod- 
ification in the sho ck jump conditions because of the presenc e 
of magne t ic fie ld s jFan et alj l2004al : IZhang & Kobavashil l2005h . 
iLvutikovl j2006l) : iGenet et all j2006l) studied the afte rglow from 
sub-fast (cr » 1) MHD models. iGiannios et all j2008l) found that 
rarefaction wav es, instead of shocks, can also form in mag netized 
ejecta (see also lMizuno et al.ll2009l : lAlov & Mimicall2008l) . These 
studies have left several issues open (e.g., that of the timescale of 
transfer of the magnetic energy into the shocked external medium) 
or debated (see ILvutikovl 120051) . The complete dynamical evolu- 
tion of homogeneo us, magnetized ejecta colliding with the ISM has 
been simulated by jMimica et al.ll2009bl MGA09 in the following) 
for a few characteristic cases, showing how the magnetic energy is 
transferred in the shocked ISM until a self-similar blast wave forms 
at later time. 



Although bolometric light curves of the expected emission of 
magnetized ejecta were computed by MGA09, a detailed calcula- 
tion of the synthetic light curves at different wavelengths is miss- 
ing. I n this work we have coupled the SPEV (SPectral £Volution) 
code jMimica et al ■l2009al) . used to compute the non- thermal emis- 
sion from re lativistic outflows, to the MRGENESIS jMimica et al.l 
I2005ll2007l) RMHD code in order to study in accurately the signa- 
tures of magnetic fields and the (non)existence of shocks in after- 
glow ejecta. 

The paper is organized as follows: in Section[2]we give a brief 
overview of the results in MGA09 and present our extended set of 
models. Section [3] gives details of the model we use to compute 
the non-thermal synchrotron emission from relativistic shocks. In 
Section [4] we discuss the energy content of GRB ejecta and of the 
external medium as a function of the distance and of the initial flow 
magnetization. Light curves and spectra are presented in Section[5] 
We apply our models to the GRBs 990123 and 090102 in Section[6] 
Discussion and conclusions are given in Section|7] 



2 NUMERICAL MODELING OF EJECTA-MEDIUM 
INTERACTION 

We are interested in following the evolution of the flow well af- 
ter the collimation, acceleration, and GRB emission phases. At this 
stage, the ejecta are moving radially, cold (because of adiabatic ex- 
pansion) and expected to contain a dominant toroidal component of 
the magnetic field. We parametrize the magnetization of the ejecta 
through o"o = B'^ c /4npc 2 , where B[ nac and p are the rest-frame 
macroscopic (ordered) magnetic field strength and density respec- 
tively. 

We focus on the initial stages of the interactions during which 
only moderate deceleration of the ejecta takes place. The ejecta are 
expected to have bulk Lorenz factor r » 1/8, where 9 is open- 
ing angle of the jet. This allows us to ignore 2D effects (such as 
lateral jet spreading) justifying our assumed spherical symmetry. 
The ejecta is modeled as a homogeneous shell of initial thickness 
Ao <k ro, where ro is the initial distance determined by the fol- 
lowing considerations. Because of its magnetic pressure, the shell 
will tend to expand developing rarefactions on both radial edges. 
These rarefactions propagate at fast magnetosonic speed in the 
ejecta. Comparing the time it take s for a rarefaction to cross the 
ejecta with the expansion timescale. lGiannios et al. I j2008l) estimate 
the distance at which the rarefactions broaden the shell. Further- 
more, a reverse shock may form in the ejecta as result of the in- 
teraction with the ISM, or as the late steepening of the conditions 
in the tail of the primordial rarefaction sweeping backwards the 
ejecta from its leading radial edge (MGA09). We set the initial dis- 
tance r of our study such that it is large enough to make the study 
amenable to numerical RMHD simulations, but small enough so as 
to begin our simulations in the coasting phase of the ejecta, before 
the deceleration of the ultra-relativistic flow. Hence, we can track 
self-consistently, both the evolution of the complete Riemann fan 
emerging from each radial edge of the shell and the initial deceler- 
ation of the ejecta. 

The single quantity that determines the strength of the re- 
verse shock in unmagnetized ejecta is £ = V'/A)F 4/3 where 
I = (3E/4mi ext m p c 2 ) 113 is the Sedov length, zi cx( is the number den- 
sity of the external medium, c is the speed of light, m p is the pro- 
ton mass, and E = Et = 47irj 2 .Aor 2 ja s heiic 2 is the total (kinetic, Et) 
energy of the ejecta (the internal energy contribution is neglected, 
since the shell is initially assumed to be cold). The £ quantity can 
be related to the more familiar "deceleration radius" rd cc defined 
as the radius where the ejecta accumulate a mass r^ 1 times their 
own mass from the external medium, resulting in r dec = //T? . 
For parameters relevant for GRB flows, the parameter £ is expected 
to be of order of unity (within one order of magnitude). The sec- 
ond parameter that naturally appears when studying the decelera- 
tion of magnetized ejecta is the magnetization ctq in the coasting 
phase (cro is much less constrained from theory or observations). 
We note that in magnetized ejecta, the total (kinetic plus magnetic) 
energy is E = Et(l + cro). Realistic GRB models involve ejecta ini- 
tially moving with Lorenz factor To ~ 100 - 1000 and with thick- 
ness A ~ cT G rb ~ 10 12 cm <k r dec (T Gm is the duration of the 
GRB phase). Such simulations pose large numerical challenges. In 
MGA09 we overcome these challenges by showing that one can 
simulate a relativistically moving shell of, say r ~ tens and pre- 
dict the evolution of a shell with larger To but with the same <xo and 
f . We refer to this procedure as rescaling. 

In this work, we extend the models of MGA09 to explore a 
wider parameter space of <t for two values of £ = 0.5 and f = 1.1. 
We refer to the £ = 0.5 runs as "thick shell" and f = 1.1 as "thin 
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Table 1. Properties of the simulated ejecta-external medium interactions. 
See text of Sect.[2]for the definitions of parameters. 

thin shell = 1.1) thick shell (£ = 0.5) 



E [erg] 


3.33 x 10 53 


3.33 x 10 54 




0, 0.0001, 0.001, 0.01, 0.05, 0.1, 0.5, 1 


«ext [cm" 3 ] 


10 




Pexl/Pext c ~ 


io- 3 




P 'shell /PshellC 2 


2.5 x 10" 


3 


r [cm] 


5 x 10 16 




'"O.rescaled [cm] 


6.8 x 10 15 


3.7 x 10 15 


A [cm] 


10 15 


10 16 


Ao,rescaled [cm] 


3.4 x 10 11 


3 x 10 11 


r 


15 




r"o,rescaled 


300 


750 



shell" models respectively. In Tab. Q] we summarize parameters of 
the numerical models. E, <x , r , A and T are the shell total energy, 
magnetization, initial radius, width and initial Lorenz factor, re- 
spectively. We assume that both the external medium and the ejecta 
are ideal gases which obey an equation of state with a fixed adia- 
batic index y a( j = 4/3, that p = nm p and give in Tab.Q]the comov- 
ing number density of the external medium « cxt , and the pressure 
to rest-mass density ratios P ext /p ex ,c 2 and /'shell /Psheiic 2 in the exter- 
nal medium and in the shell, respectively. Finally, we indicate the 
Lorenz factor ro^escaied and shell thickness Annealed to which results 
of our simulations are rescaled (see Section[5]l. 

For the numerical simulations we use the relativistic mag- 



netohydrodynamic code MRGENESIS dMimica et ai1l2005l. 12007 
a high-resolution shock capturing sch eme based on GENESIS 
dAlov et alJl 19991 : iLeismann et alj|2005t) . In the problem at hand, 
we assume spherical symmetry and the fluid is discretized in spher- 
ical shells (zones). As described in detail in MGA09, we use the 
grid re-mapping technique, which enables us to decrease the com- 
putational cost of our simulations by concentrating the computation 
on a domain, which is several initial shell widths wide, and which 
follows the shell as it moves radially outwards. Since the disadvan- 
tage of this is that at later times we lose a part of the tail of the 
ejecta, we have recomputed some of the models listed in Tab. [T] 
with a five times wider radial grid. Thus, we can correctly follow 
the evolution of the total magnetic energy content to later times (see 
Section|4ll. 



3 NON-THERMAL EMISSION 

The bulk of the observed radiation that follows the prompt GRB 
phase is believed to be synchrotron emission from electrons ac- 
celerated at the FS and at the RS. Behind the shocks, we as- 
sume that particles are accelerated to a power law distribution with 
N(y) oc y~ p ', where p is the power-law index of the electron distri- 
butiorQ. 



1 Recent particle-in-shell (PIC) simulations iSpitkovskvl 1200 8|) in- 
dicate that a pronounced Maxwellian electron component appears 
around the minimum of the non-thermal distribution in unmagnetized 
shocks. Observational c onseq uences of this component are discussed in 
iGiannio s & Spitkovsky 12009). Furthermore, strong, large scale magnetic 
fields inhibit the formation of a no n-thermal tail in the particle distribution, 
as found by the PIC simulations of lSironi"& Spitkovsky (2009j). In Section 



To compute the minimum Lorenz factor ymin of the electron 
distribution, we assume that a fraction e c of the energy dissipated 
by the shock is used to accelerate electrons. We use the same e e 
parameter for both the RS and the FS. We estimate the maximum 
Lorenz factor y max to which electrons are accelerated by equating 
the synchrotron cooling time with the gyration time. Then we have 



_ 3P sh mp p-2 

Tmin t « Q 



p sh c A m e p ■ 



1 



a 2 4 \ '/2 
Am r c ' 



47re 3 fi' 



(1) 



where P^ and p s i, are the thermal pressure and comoving density 
of the fluid downstream of the shock, m e is the electron mass, and 
e is the electron charge. 

Shock compression of the typical ordered external medium 
magnetic fields is not s ufficient to yield a subs tantial synchrotron 
emission (see, however, iKumar & Duran|[2009l and Sect. 6 of this 
work for possible exceptions) and a stochastic field B' st amplified 
by plasma instabilities in the shock is assumed. We model such a 
stochastic field in a standard way, by assuming that its energy is a 
fraction e B of the internal energy of the fluid behind the FS, i.e. the 
comoving stochastic magnetic field strength is B' sl = -\j24n€ B P sh . 
The same instabilities may operate in the reverse shock (if the ini- 
tial magnetic field of the ejecta is weak). On the other hand, in high 
<t models, the shock compressed field dominates. Thereby, in or- 
der to have a continuous description of the emission of magnetized 
and non-magnetized models, for the purpose of computing the syn- 
chrotron emissivity, we consider that the total comoving magnetic 
field strength is the sum of the stochastic field plus the ordered 
macroscopic field, i.e., B' = B' m , x + B' st . 

Using the SPEV code, we calculate the synchrotron emission 
from the injected particles, and also follow the evolution of the par- 
ticle distribution downstream the shocks including cooling (adia- 
batic a nd radiative) and synchrotron self-absorption dMimica et al.l 
l2009al) . Our calculations neglect synchrotron-self-Compton (SSC), 
whic h may make sign ificant contributions in the high-energy emis- 
sion dFan et all l2008h and, in some cases, it can also affect the 
strength of the synchrotron emission because of intensive Com p- 
ton cooling of the electrons dPetropoulou & Mastichiadisll2009h . 

We note that SPEV tak es a different a pproac h from 
other methods such as that of Ivan Eerten et alj d20ld) . While 
Ivan Eerten et alj d2010h perform a part of the emission calcula- 
tion during the RMHD simulation, SPEV operates on the post- 
processed data of the RMHD simulation. To model the spa- 
tial distribution of non-ther mal particles it uses a volume-filling 
Lagrangian particle scheme dMimica etal]|2009al Appendix A). 
We modified SPEV so that, if the synchrotron absorption is ne- 
glect ed, it compu t es the flux density according to the equatio n 
3 in iGranot etal] dl999i) (see also equation 4 in iDermerl [20081) . 
In case when synchrotron absorption is not neglected, we solve 
the radiative transfe r equation by following photon trajectories 
dMimica et alj|2009al . Appendix A). 



4 DYNAMICS AND TIME SCALES 

The dynamics of the deceleration of the ejecta has been described 
in detail in MGA09 for <x = and cr = 1. Here, we have run 
several additional simulations to investigate intermediate values of 
ctq. Since the methodology to compute the emission is much more 



5, we comment on the observational implications of the lack of particle ac- 
celeration. 
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elaborated in the present paper (Sect. 0, we have validated and 
refined the MGA09 findings (regarding the observational signature 
of the models). The main results of MGA09 can be summarized as 
follows: 

• The deceleration depends on just two parameters £ and ctq. 
Lower £ and <x result in stronger dissipation of energy in the re- 
verse shock and vice- versa. 

• For the typical parameters of the GRB flows, cr ~ 1 strongly 
suppresses the amount of energy dissipated in the reverse shock. 
The front part of the ejecta initially forms a reverse rarefaction. At 
a later stage, because of the spherical expansion, a weak reverse 
shock forms at the rarefaction tail. 

• The total magnetic energy of the ejecta remains approximately 
constant or even increa ses while they are crossed by the reverse 
shock (also discussed in lZhang & Kob avashi 20051). 

• After the RS reaches the rear end of the ejecta, a rarefac- 
tion forms and crosses the ejecta and the shocked external medium 
within the expansion timescale of the blast wave. The rarefaction 
dilutes the magnetic energy in the ejecta. When the rarefaction head 
catches up with the FS there is a moderate but sharp decrease of 
the Lorenz factor of the shock. Only at this stage the solution ap - 
proaches the self similar evolution of lBlandford & McKed dl976t) . 

A common feature of all our simulations is that the self-similar 
solution is established at an observer time that is longer than the du- 
ration of the burst (which corresponds to the light-crossing time of 
the shell thickness). At earlier times, there is substantial energy in- 
jection at the FS that clearly affects the rate of deceleration of such 
shock (see Sect.|5]for lightcurves and Sect. 7 for observational im- 
plications of this effect). The self-similar evolution can, therefore, 
be assumed only when calculating afterglow lightcurves well after 
the duration of the burst. 

While we cannot exclude that there exists a region of the 
- (To) parameter space in which the energy transfer takes place 
before the end of the GRB, we consider this possibility unlikely. 
The reasoning is the following. It can be shown that the observer 
time corresponding to the reverse shock cr ossing of the ejecta i s 
larger than the duration of the burst (e.g., iNakar & Piranll2004h ; 
magnetization o f the ejecta does not change this conclusion (e.g., 
iFan et alj|2004ah . At this stage a rarefaction forms and crosses the 
shocked ejecta. Our simulations show that energy is transferred at 
the FS until the rarefaction reaches the FS. Since the rarefaction 
moves slower than light, the energy transfer is completed after the 
peak of the reverse shock emission (which trails the bulk of the 
GRB emission). 



4.1 The long-term evolution of the magnetic energy 

Although, as we just discussed, the bulk of the energy of the ejecta 
is passed onto the shocked ISM on a short time scale (slightly 
longer than the GRB duration), there is a residual magnetic energy 
remaining in the ejecta. Such energy may, however, have important 
observational i mplications if di ssipated at later times. For example, 
it is argued in iGiannio si d2006h that the slowing down of magne- 
tized ejecta leads to revival of (current driven) MHD instabilities. 
The instabilities can lead to the dissipation of magnetic energy in 
localized reconnection regions which may pow er the observed af- 
terglow X-ray flaring ^Burrows & et alj|2005l) . In this picture no 
late central engine activity is needed. Our ID models cannot tackle 
the issue of the nature of such (3D) instabilities, but they can quan- 
tify how much energy remains in the magnetized shell as function 
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Figure 1. Magnetic energy fraction as function of the dimensionless ob- 
server time tco ■ Thin lines show results for thin shell models o"o = 0.1 and 
o"o = 1 with the short grid, while the thick lines show results for the long 
grid. The rapid drop of the magnetic energy at late times is not physical but 
due to the fact that the whole magnetized shell leaves our computational 
box. This problem is ameliorated in the long grid models, where it shows 
up at later times, when the magnetic energy content is much smaller. 



of observer's time that could potentially power any afterglow flar- 
ing. 

In Fig.[T] we plot the magnetic energy as function of a dimen- 
sionless observer time defined as t CD = (ct- R C d)I where R C d is 
the location of the contact discontinuity separating the shocked cir- 
cumburst medium from the ejecta. Since any emission from mag- 
netic dissipation takes place behind the contact, / C d gives the short- 
est time after the burst that any emission from the ejecta will be 
observed. 

The remapping technique we are using to track the evolution 
of the shell consists in adding a certain number of numerical zones 
in front of the leading forward shock. The same number of zones is 
removed from the rear part of the domain. Since the shell is ultra- 
relativistic, the rear boundary is causally disconnected from the RS 
during the time in which we evolve our models. Thus, even consid- 
ering that a part of the shell is lost through the rear radial bound- 
ary, for the purpose of computing the dynamics (or the emission) 
of the leading forward and reverse shocks, it suffices a relatively 
small computational domain L spanning not only the shell, but a 
small portion of the leading and of the trailing medium, namely 
L = 1 .5Ao (L = 4Ao) for the thick (thin) shell models. However, due 
to the loss of the rear part of the shell, the calculation of any integral 
quantity extending over the whole shell (in particular the magnetic 
energy) is underestimated. To ameliorate such a deficiency of our 
models, we have computed some of the "thin shell" simulations 
with a ~ 3 times larger grid ("long grid" in Fig. [TJ than that used 
in MGA09 ("short grid" in Fig.[TJ. At ?cd < 1 the magnetic energy 
changes little (and actually increases in the <x = 0. 1 run because 
of the (reverse) shock compressing the ejecta). The magnetic en- 
ergy drops fast at tco ~ 1 because of the rarefaction that follows. 
Nevertheless, a fraction of ~ 10% of the total energy remains in 
magnetic form for longer time showing only modest decline with 
time (note that the large grid is essential to follow this late evolu- 
tion). We find that at least ~ 2% of the energy in still in magnetic 
form at t CD ~ 10 (i.e., about 10 times the burst duration). The rapid 
drop of the magnetic energy at late times is not physical. It is due to 
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Figure 2. R-band light curve for the thin shell with initial magnetization 
ctq = 0.01. Shown are the light curves originating from the FS (thin full 
line), RS (dotted line) and the combined light curve (thick full line). In the 
calculation eg = 5 X 10~ 3 , e e = 10 - ' and the electron power-law index 
p = 2.25. The shaded area denotes the observational interval which is si- 
multaneous with the prompt emission. 

the fact that the contact discontinuity reaches the back of our grid 
and the whole magnetized shell leaves our computational box (and 
we are left with the shocked ISM). 

In view of these results, we suggest that there is sufficient en- 
ergy remaining in the magnetic field for at least ~ 10 times the 
duration of the burst to yield a substantial emission if a large part 
of the magnetic energy is dissipated and radiated away. 



5 LIGHT CURVES AND SPECTRA 

In this section we present the non-thermal light curves and spectra 
computed from our models. We first discuss the generic light curves 
corresponding to thin and thick shells for different magnetization of 
the ejecta. We also show a prototype emission spectrum at the end 
of this section. 

5.1 Thin shells 

Figure [2] shows the R-band light curves for the thin shell (i.e., £ = 
1.1) model with initial magnetization tr = 0.01 and bulk Lorenz 
To = 300. The emission from the shocked external medium (thin 
full line), shocked ejecta (dotted line) and the total emission (thick 
full line) is displayed. We have assumed e B = 5 x 10~ 3 and e e = 
10" 1 so as to calculate emission. The total R-band light curves for 
models with different values cr from to 1 are shown in Fig. [3] 

As can be seen in Fig. [2] and, in general, in Fig. [3] for the pa- 
rameters chosen the RS emission is marginally brighter than that 
of the FS in optical bands. Moderate differences in light curves for 
t<25 s are caused by the effect of <t on the RS emission. Increasing 
(To from to 0.01 leads to progressively brighter RS emission be- 
cause of increased synchrotron emissivity resulting from a stronger 
field. Further increase of the magnetization greatly weakens the RS 
emission and the amount of heating that takes place at the shock 
(and, therefore, reduces the electron minimum Lorenz factor in the 
downstream fluid; see Eq.QJ resulting in a negligible optical emis- 
sion compared to that of the FS. 
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Figure 3. Combined (FS + RS) R-band light curves for thin shell models 
with (T = 0, 10~ 3 , 10~ 2 , 10~' and 1 (thin full, dotted, thick full, thin 
dashed and thick dashed lines, respectively). The rest of the parameters are 
like in Fig. [2] 
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Figure 4. Forward shock emission. R-band, 1 keV and 100 MeV light 
curves (full, dotted and dashed lines, respectively) for thin shell models 
with o"o = (thin lines) and <tq = 1 (thick lines). 

The FS emission shows a weak dependence on <ro, the differ- 
ences being more pronounced at higher energies (X-rays, gamma- 
rays) as shown in Fig. [4] Our models display an achromatic break 
in the light curves shortly after the end of the burst. The break takes 
place at t ~ 30 s and 20 s for the <xo = and (To = 1 models, re- 
spectively, and it results from the rarefaction that reaches the FS and 
decelerates it. At early times, the emission in gamma- and X-rays 
is systematically brighter and flatter for high-cro models than in the 
optical band, while the gamma- and X-ray flux rises faster for low- 
<To models. Observing early enough after the burst, the high-energy 
emission from the FS can contain a valuable information about the 
magnetization of the ejecta. If the afterglow is observed too late, 
in this example after approximately 30 s(Fig.|4j. light curves with 
radically different values of a~o become indistinguishable, since the 
ejecta enters into the self-similar evolution phase, dependent only 
on the total initial energy and the external medium density, which 
are the same for all models. 
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Figure 5. Same as Fig. [2]but for the thick shell with initial magnetization 
<r = 0.01. 

5.2 Thick shells 

We, now turn to the emission from a thick shell = 0.5) for To = 
750. Figure [5] shows the contributions from the FS and the RS to 
the total R-band light curve for the model with initial magnetization 
(To = 0.01. We can see that the RS the contribution to the emitted 
flux here is much more important than in the thin shell case. The 
RS dominates the emission until t = 50 s. 

As can be seen in Fig. [6] the RS emission dominates in the 
early afterglow over that of the FS in the optical bands. Large dif- 
ferences in light curves for t < 50 s are caused by the effect of 
the magnetization <x on the RS emission. For <x = 0, there is a 
clear RS signature. It results from the synchrotron emission of non- 
thermal particles gyrating in the stochastic magnetic field (presum- 
ably) generated at the shocks. Increasing cr from to 0.01 - 0.1 
leads to brighter RS (synchrotron) emission because of the stronger 
(shock compressed) magnetic field. Maximal RS emission happens 
for cr ~ 0.01 - 0.1. Around such values of <x there is a turnover 
of the emission, which drops as the magnetization keeps growing. 
In effective terms, the RS light curve changes only moderately for 
10~ 3 < cr < 0.1, and hence, inferring an accurate value of <xo from 
a careful modeling of the RS emission is difficult. Further increase 
of the magnetization greatly weakens the reverse shock and sup- 
presses the energy injected to the particles. As a result, there is 
a clear suppression of the reverse shock emission for <r ~ 1 or 
higher. 

Other mild differences in the optical lightcurves are caused by 
differences in the Lorenz factor of the FS at early times because of 
the different <x as shown in Fig. [7] The thick-shell models show 
the same achromatic break in the light curves shortly after the end 
of the burst as seen in the thin shell ones. In this example, after ap- 
proximately 30 s, FS light curves become indistinguishable, since 
the ejecta evolution approaches to the self-similar evolution phase. 

5.3 Emission spectra 

The instantaneous spectrum of the thick shell model <x = 0.01 
taken at / = 10 s is shown in Fig. [8] This time approximately corre- 
sponds to the peak of the RS emission and to the end of the prompt 
(~MeV) emission. Thick full line shows the total spectrum, while 
thin full and dashed lines show the FS and RS contributions, respec- 
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Figure 6. Same as Fig. [3] but for thick shell models. 
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Figure 7. Same as Fig. [4] but for thick shell models. 

tively. Both spectra have been computed neglecting the synchrotron 
self-absorption (SSA). Thin dotted line shows the RS spectrum 
when SSA is taken into accounj^]. 

The peak of the vL v emission produced by each of the shocks 
corresponds to their respective cooling break frequencies. It ap- 
pears in the X-ray band and in the UV band for the FS and RS, 
respectively. In Fig.[8]we also show the frequency of the SSA break 
for the reverse shock emission (v a ,Rs) that lies in the infra-red. Not 
surprisingly, the optical emission in the range ~ 10 14 - 10 15 Hz is 
dominated by the RS, while the X-ray-to-y-ray emission is domi- 
nated by the FS. 

In deriving these results, we have assumed that electrons are 
accelerated into a power law distribution at the shocks. Despite the 

2 We do not include SSA in all the light curve calculations since the com- 
putational costs of such a calculation are prohibitive (one of the main rea- 
sons is that the algorithm which includes SSA is not as easily parallelizable 
as the optically-thin calculation). However, we have verified that the light 
curves of all the bands shown in this work lie above the SSA frequency 
v a ,RS, making the SSA correction negligible for the purpose of the present 
paper. 
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Figure 8. Emission spectrum at t = 10 s for the thick shell model with 
o"o = 0.01. Thin, full and dashed lines show optically thin FS and RS spec- 
tra, respectively. Thick full line shows the total spectrum. Thin dotted line 
shows the portion of the RS spectrum which is optically thick to synchrotron 
self-absorption (see text for details). Vertical dotted lines denote FS and RS 
cooling break frequencies v c fs and v c rs , as well as the RS synchrotron-self 
absorption break frequency v a ,Rs ■ 



fact that particle acceleration in (relativistic) shocks is poorly un- 
derstood, there are indications that strong, large scale fields along 
the shock plane can inhi bit particle acceleration. Fo r instance re- 
cent PIC simulations by ISironi & Spitkovskvl J200gh have shown 
that for <x > 0.01, magnetized, relativistic shocks result in efficient 
heating of the particles with the distribution lacking non-thermal 
tails. If this result holds for the mildly relativistic shock forming 
in the ejecta, such lack on high-energy tail can substantially af- 
fect the RS emission (and, as a matter of fact, any emission from 
internal shocks). If the characteristic frequency (associated to the 
thermal distribution of the particles) lies bellow the observed band 
then there will be practically no observed emission from the RS. 
When the characteristic frequency of the RS emission (as is often 
the case) lies in the infrared, bursts of sufficiently high cr ^0.01 will 
be optically dim (there are no electrons accelerated to high enough 
energy to (synchrotron) emit in the optical or harder bands). This is 
an additional potential mechanism that suppresses the RS emission 
in moderately to strongly magnetized ejecta. 



6 APPLICATION TO GRB 990123 AND GRB 090102 

In this section we apply our early afterglow models to two specific 
GRBs exhibiting bright optical emission which may be powered by 
a reve rse shock, GRB 990 123 dBriggs et al|[l999llAkerlof & et all 
1 19991) and GRB 090102 dGendre et al.l 120091 : ISteele et alj|2009l) . 
The famous GRB 990123 had a 9th magnitude optical flare a few 
seconds after the end of the GRB and is one of the brightest ever 
detected. 090102 is a rather bright GRB wh ose early optical af - 
terglow shows declining, but bright emission dGendre et ai1l2009l) 
and, most importantly, a significant ~ 1 0% polarization, in dicating 
large scale fields in the emitting region dSteele et al.l2009l) . 

Several of the parameters of our model can be directly inferred 
by the observed properties of these bursts. The rest-frame duration 
of both bursts is T GRB as 10 s, constraining the shell thickness to be 
A - cT GRB = 3x 10 11 cm. The gamma-ray fluence of GRB 990123 



and 090102 is E y =s 3 x 10 54 erg and E 7 - 6x 10 53 erg, respectively. 
Assuming ~ 20% efficiency for the prompt emission mechanism 
the total (isotropic equivalent) energy of the ejecta can be estimated 
to be £990123 - 1.5 x 10 55 erg and £090102 - 3x 10 54 erg. Finally, the 
fact the the reverse shock peaks very shortly after the burst indicates 
that £ < 1 in these bursts. 

Taking these observational constraints into account, we look 
for the values To and <To for which our rescaled numerical models 
best fit the observations. Figure|9]shows the R-band light curves of 
GRB 990123 (circles) and GRB 090102 (trianglesfl We also plot 
the light curves of two thick shell models, the <xo = 0.01 rescaled 
to T = 640 (thick lines) and cr = 0. 1 rescaled to T = 940 (thin 
lines). Full, dot-dashed and dashed lines denote the total, RS, and 
FS emission. The data in our simulation only allowed us to compute 
the FS emission up to t = 80 s. By this time, however, the evolution 
of the blast is self-similar. For the purpose of better visualization, 
we have used the Blandford-McKee asymptotic solution to extend 
the FS light curve to later times. 

Even though our best fits suggest that ctq = 0. 1 and <tq = 0.01 
for GRBs 090102 and 990123 respectively, <r is not well con- 
strained. From Fig. [6] one can see that the RS emission varies rela- 
tively little for a large range of 0.01 < <xo < 0. 1 indicating that pin- 
pointing the magnetization with precision from fitting the reverse 
shock light curve is hard. 

The optical afterglow of GRB 090102 shows ~ 10% polar- 
ization at an exposure th at took place 160 - 220 s after the GRB 
trigger dSteele et alj2009h . This value is very hard to understand as 
synchrotron emission from small scale magnetic fields but is still 
smaller than t hat expected from a coherent field. The exposure time 
interval of the ISteele et all j2009h observations corresponds to the 
time interval 100 - 140 s in the frame of the source. One can see in 
Fig.|9]that at ~ 100 s the FS emission may already dominate over 
the RS emission by a modest factor. The observed 10% polarization 
may, therefore, be understood as a combination of the very weekly 
polarized, brighter FS emission and a strongly polarized emission 
coming from RS that propagates in a flow of coherent field. If this 
interpretation is correct, earlier polarimetry of afterglow emission 
should reveal even stronger polarization signal, for these few bursts 
for which the RS emission is sufficiently bright with respect to the 
FS. 

Another intriguing result from the fitting of our models is that, 
in order for the FS to account for the late-time emission, a rather 
small value of e B ~ 10~ 7 needs to be assumed. For GRB 9 90123, 
a similar finding is presented in lPanaitescu & KumaJ d2004l) . Note 
that such weak fields may be simply due to shock compression of 
the ISM magnetic fields (i.e., there is no need for field generation at 
the shocks for these bursts). I Kumar & Duranl d2009l) present more 
bursts whose afterglow modeling is compatible with very low val- 
ues of e B - 



7 DISCUSSION AND CONCLUSION 

The theoretical modeling of GRB central engines is still far from 
predicting whether the GRB flow is thermally or magnetically dom- 
inated in the launching region. Furthermore, prompt GRB emission 
mechanisms are not understood well enough to enable probing of 
the magnetization of the emitting plasma. A more promising way 

3 To convert flux to luminosity, we have assumed for both GRBs a redshift 
Z = 1.6, and a standard cosmology model with the Hubble and density 
parameters are Hq = 71 kms~' Mpc~', Cl,„ = 0.27, and Cl\ = 0.73. 
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Figure 9. R-band li ght curve for G RB 990123 iBriggs etai] Il999l; 
Akerlof & et al.1 1 19991) . GRB 090102 jGendre et alj 1200ft ISteele et all 
20091) jlotted in the rest frame of the central engine of the burst. Shown 
are the thick shell models which most closely match the observed light 
curves, assuming for both bursts ejecta width A„ = 3.3 X 10" cm (with 
corresponding duration of 1 1 seconds; marked by the shaded region in the 
plot). The ctq = 0.01 model (thick lines) has been rescaled to Tq = 640 
and n cxt = 10 cm -3 . The o"o =0.1 model (thin lines) has been rescaled to 
T = 940, ;i cxt = 1 cm 4 . For both curves e e = 0.02 and e B = 4 X 10~ 7 
are assumed. Full, dot-dashed and dashed lines denote the total (FS + RS), 
RS and, FS emission. The vertical dotted line shows the observer time up to 
which the data from the simulation has been used to compute the FS emis- 
sion; from that point onwards the Blandford-McKee asymptotic solution 
has been used to compute the FS light curve. 



to infer the magnetization of the GRB flow may come through 
studying interactions with the circumburst medium. Quantitative 
inferences, however, critically depend upon the detailed model- 
ing of the dynamics of such interaction, as well as the emission 
mechanisms from shocked plasma. In this work we presented rel- 
ativistic MHD models of ejecta-medium interactions coupled to a 
radiative transfer code to calculate multi-wavelength synchrotron 
lightcurves. The most salient results of our parametric study are the 
following: 

(i) For very weakly ma gnetized ejecta, t he standard for- 
ward/reverse shock structure ( I Sari & P iran 1995) is reproduced by 
our simulations. The reverse shock can dominate the early opti- 
cal/IR emission while the forward shock dominates the subsequent 
emission at all wavelengths. Increasing the magnetization of the 
ejecta to cr ~ 0.01 -0.1, the reverse shock remains strong while 
the powerful fields result in a brighter synchrotron emission, typ- 
ically, in the optical band. At still stronger magnetization cr > 1, 
the dissipation at the reverse shock and, therefore its emission are 
greatly reduced. 

(ii) Observationally, the maj ority of GRBs do n ot show bright 
optical emission at early times. iKlotz et ail d2009h present a sys- 
tematic study of early optical observations and tight upper limits 
tens of seconds after GRB emission. We interpret this lack of early 
optical detections as evidence for weak/absent reverse shock as ex- 
pected from GRB flows characterized by cr > 1 at large distance 
from the central engine. Nevertheless, a few GRBs such as 990123 
and 090102 do show bri ght optical emission that is likely coming 
from the reverse shock I Akerlof & et al .11 19991 ; iMeszaros & Reesl 
ll999l ; lNakar & Piranll2005l ; ISteele et al.ll2009l) . 



(iii) Our detailed modeling of early afterglow emission from 
GRB 990123 and 090102 indicates that the magnetization of the 
ejecta is <r ~ 0.01 -0.1 for these bursts, i.e., close to the value that 
maximizes the RS brightness (see Fig. [9>. For GRB 090102, the 
presence of a large scale field is s trongly supported by the detec- 
tion of ~ 10% optical polarization dSteele et alj|2009t) during a 60 s 
exposure shortly after the burst. It, therefore, appears that there is a 
range of magnetizations of tr > 0.01-0. 1 for the GRB ejecta at the 
radius of deceleration. This values are compatible with super-fast 
magnetosonic MHD models for GRBs. 

(iv) The reverse shock slows down and compresses the shell. 
As a result, during the RS crossing the magnetic ene rgy of the 
ejecta is actually increased (Zhang & Kobayashi 2005, MGA09). 
Subsequently, the bulk of the magnetic energy is diluted because 
of a rarefaction that crosses the ejecta (produced when the RS 
reaches the innermost radial edge of the ejected shell). The remain- 
ing (~10%) of the magnetic energy decays slowly with time (see 
Fig. 1). If a substantial fraction of this energy is radia ted away be- 
cause of dissipative MHD processes jGiannio it may ac- 
count for the X-ray flares ob served minutes to hours after many 
GRBs l lBurrows & etai]|2005h . 

(v) The bulk of the magnetic energy is transferred into the 
shocked ISM shortly after the end of the burst. Only at this point, 
the blast wave establishes the well known BM self-similar behavior. 
A flat or rising light curve until well after the GRB is over at high- 
energies is the signature imprinted on the FS emission by the afore- 
mentioned energy transfer. An initially flat high-energy light curve 
may be the signature of <x ~ 1 ejecta while a rising one comes 
from smaller magnetization (assuming always a constant density 
ISM). When the energy transfer is completed, an achromatic break 
takes place. After the break, the light curve declines as expected 
from the BM solution. 

(vi) It has bee n proposed jKumar & Durai] 120091 : 
iGhisellini etai] l20ld) that the GeV emission observed by the 
LAT detector on FERMI in several bursts is result of the FS. The 
model can account for GeV observations shortly after the end 
of the burst and and for the late optical and X-ray data. In its 
simplest form, however, the FS model cannot explain the peak of 
the GeV emission at times t peak < r G RB seen in GRB 080916C, 
090510 and likely in 090902B (since, as we have shown here, the 
energy tr ansfer to the FS r esults in t pea k > ?grb)- On the other hand, 
from the IGhisellini et alj feOld) sample of several LAT bursts, at 
least three GRBs (090323, 090626, and 091003) show a rather 
fiat GeV light curve that steepens only after the MeV emission 
(i.e. the traditional GRB) is over. Though more detailed modeling 
is needed, the emission from these bursts appears compatible 
with coming from the early FS driven by ejecta with cro > 1 in 
agreement with inferences we have made from the lack of observed 
of RS emission. 



The results outlined so far suggest that to explain the paucity 
of observed optical flashes one needs to argue that most of the ob- 
served bursts posses substantial magnetizations (point (ii) above). 
However, our models are sensitive to a number of parameters, 
among which, e B is probably the most crucial one. The value that 
we use as representative in our models (e B ~ 5 x 10~ 3 ; Sect. [5} re- 
sults in powerful optical flashes even for the cr = model. Values 
of e s <k 5 x 10~ 3 reduce the strength of the RS emission. We find 
that for e B <IQ 6 - 10~ 5 the RS emission falls below that of the FS, 
for all our purely hydrodynamic models (cr = 0). Therefore, even 
in fireball models one may argue that the optical emission of the RS 
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is suppressed if mildly relativistic shocks are not able to efficiently 
generate stochastic magnetic fields. 

Another issue that has to be taken into consideration is the fact 
that we are modeling the ejecta as initially being cold. Moderately 
hot or warm ejecta (i.e., in which Pipe 1 < 1) may result, e.g., if the 
conversion of internal to kinetic energy is slower than predicted by 
the standard fireball model (which might not be unlikely by look- 
ing at the numerical models of, e.g.. lAlov et al.|[200ol 12005^. or if 
internal shocks heat up the ejecta shortly before the reverse shock 
crossing. In such scenario, a weaker RS may be envisioned, since 
it is more difficult to shock warm than cold matter. Thus, we may 
not exclude the possibility that the absence of a clear RS signature 
in many GRBs arises (assuming that the ejecta is not magnetized) 
from the fact that the ejecta is warm by the time the afterglow emis- 
sion takes place. 
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